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Abstract 

The eigenstates of a single massless Dirac particle confined in a linear potential are calculated 
exactly by direct solution of the Dirac equation. The electromagnetic structure functions are 
calculated from the Dirac wave functions of the ground and excited states of the particle by coupling 
to its conserved vector current. We obtain the longitudinal and transverse structure functions as 
a function of y = v — |q|, where v and |q| are the energy and momentum transferred to the target 
in its rest frame. At values of |q| > 2.5 GeV, much larger than the characteristic energy scale ~ 
440 MeV of the confining potential, the response exhibits y scaling, a generalization of Bjorken 
scaling. We compare the exact structure functions with those obtained from the ground state wave 
functions in the plane wave impulse approximation. The deviation from the Callan-Gross relation 
is compared with the parton model prediction. 

PACS numbers: 13.60.Hb,12.39.Ki,12.39.Pn 
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I. INTRODUCTION 



The present work is motivated by the study of the effects of interactions on the response 
of systems of confined relativistic particles. Here we consider a single massless Dirac par- 
ticle confined by a linear potential. The Dirac Hamiltonian is solved exactly for a large 
number of eigenstates and the electromagnetic structure functions are calculated by cou- 
pling a charged leptonic probe (such as the electron) to the conserved vector current of 
the target, ^ q (x)j fJi i/;q(x). This model neglects some of the important physics such as the 
contribution of qq pairs from the vacuum and the effects of radiative gluon corrections. It is 
nevertheless useful to study the role of interactions on valence quark structure functions. In 
a future paper we will use this model to investigate the role of interactions in determining 
spin-dependent observables. l| 

The model may be viewed as a heavy-light meson, such as iu, or a baryon with a light 
quark bound to an infinitely massive diquark. In the heavy-light meson, the antiquark i and 
light quark u have opposite color charge and confinement is due to the flux tube connecting 
them. In the limit of very large mass of the i the Dirac equation for the motion of the light 
quark is exact. We neglect the weak interaction. 

Benhar, Pandharipande, and Sick 2] have shown that the transverse structure function 
of the proton Wf (y, |q|), measured in inclusive deep inelastic scattering (DIS) experiments, 
exhibits scaling in the variable y = v — |q| over a large range, 5 GeV < |q| < 100 GeV. Here 
v and |q| are the energy and momentum transferred to the target in its rest frame. The 
model described in the present work also exhibits y scaling. The y variable is related to the 
dimensionless Nachtmann variable |3| £ = —y/M, where M is the target mass. The variable 
y is appropriate to systems where the constituents of the target move relativistically. jj] It 
and the Nachtmann £ are generalizations of the Bjorken scaling variable xb- Their 

relation is given by: 

f = -V- = ^ (1) 



with Q 2 — _lq| 2 — v 1 giving £ = xb in the limit Q 2 — > oo. 

In Ref. |5| we calculated the exact response of a single relativistic scalar "quark" confined 
by a linear potential. The response of this simple model exhibited a rich behavior and allowed 
the study of a diverse range of phenomena. We observed the y scaling of the response in the 
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spacelike and timelike regions of momentum transfer, an approach to scaling qualitatively 
similar to the one found in recent inclusive experiments of DIS of electrons by proton [6J , and 
Bloom-Gilman duality [7]. We also studied various approximations to the exact result. In 
particular the on-shell approximation (OSA), which assumes that constituents of a hadronic 
target are on the mass-shell before and after interaction with the probe, admits response 
only in the spacelike region of momentum transfer due to the inequality 



v 



y^k + q| 2 + m 2 - v/|k| 2 + m 2 < |q|, (2) 



where k is the momentum of the struck constituent and m is its mass. In contrast, the plane 
wave impulse approximation (PWIA), often used in nuclear physics, treats the initial state 
constituents as bound particles, which are therefore not on their mass-shell, and Eq.(J2J) is 
not satisfied. Thus there is response in the timelike region in PWIA due to the bound nature 
of the constituents. Final state interactions (FSI), however, are neglected in PWIA and the 
state of the particle after it is struck by the probe is described as a plane wave. 

In Ref.Q we investigated the role of FSI effects on the response. There we found that FSI 
broadens the PWIA response, shifting more strength into the timelike region of momentum 
transfer. This broadening of the response due to FSI persists as |q| — > oo. The exact 
response was shown to be calculable by convolution of the PWIA response with a "folding 
function." This function describes the distribution in energy of a plane wave state in the 
linear confining well. It becomes independent of the three-momentum |q| transferred to the 
target in the limit of large |q| (i.e. the "scaling" limit). Thus the y scaling, and equivalently 
Nachtmann £ scaling or Bjorken scaling of the exact response is not a signal that FSI may be 
neglected. In fact, in the semi-relativistic model of Ref. jj] it is crucial to take into account 
FSI in order to reproduce the response even in the scaling limit, |q| — > oo. This situation 
is distinct from the case of a confining potential with a non-relativistic kinetic energy. In 
such a model, the width of the response scales linearly with |q|. This precludes the need 
to take FSI into account in the scaling limit. It is also interesting to note that there is no 
Bloom-Gilman duality in the non-relativistic model. Since the width of the response scales 
linearly with |q| the resonant structure of the response at low-|q| does not average around 
the more broad response at high-|q|. 

The semi-relativistic model leaves much to be desired. Among its shortcomings are an 
ill-defined current operator, neglect of the contribution of qq pairs, and neglect of radiative 
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gluon effects. In this work we address the most egregious of these shortcomings by working 
with the conserved electromagnetic current of the confined particles, ip q (x) r y fJt ,il> q (x). 
We proceed by solving the Dirac equation for an assumed potential 

V = V(r)hl + 0) = V{r) I ' ' | (3) 




2 

V(r) = Var. (4) 



This potential has been investigated by Page, Ginocchio, and Goldman |a| and shown to 
admit a spin symmetry which suppresses spin-orbit coupling in the hadron spectrum. This 
symmetry results in the degeneracy of states with the same orbital angular momentum £ of 
the upper component of the Dirac wave function. The states with j — £ ± -= are degenerate 
and they have Dirac quantum numbers k = — (£ + 1) and £, respectively. Recall that 
K — ±(j + |) where the upper (lower) sign is used when j = £ =f |. This potential is easy 
to work with since the lower components of the wave function are not coupled by the Dirac 
matrix of Eq. (JHJ) . 

In the next section we will briefly detail the solution of the single particle Dirac equation 
for this potential. Section lTTTl gives the calculation of the electromagnetic structure functions. 
We discuss the interpretation of the contribution to the structure functions of the negative 
energy states. Section ITVl briefly compares the exact structure functions to those obtained in 
the PWIA. In the concluding Sec. EJ we study the deviation from the Callan-Gross relation. 



II. DIAGONALIZATION IN THE STATIC BAG MODEL BASIS 

We wish to solve the Dirac equation 

H^ p = let ■ p + V] $> p = E p ^ p (5) 

nam L 1 J n/tm nn nnm \ I 

for a massless fermion in the potential given in Eq.(j2J). Here n, k, and m are the quantum 
numbers of a Dirac particle in a spherically symmetric potential. The Eq. (jSj) with potential 
given by Eqs. (|3|4|) has a discrete spectrum of bound positive energy P = + states with wave 
functions that vanish as r -> oo. However, due to the form of the Dirac matrix in Eq.(jHJ) 
the wave functions of negative energy states are unbound. 

We expand the eigenstates ^ p Km in a complete basis of static bag model (herein "bag 
model") states (The explicit form of the bag model states & p Km is given in the 
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Appendix |XJ) The expansion is 

q p = \A + $ + + A~ $- , 1. (6) 

nnm / j L una aum ' una. a(— K)mi V / 



a=l,...,oo 



Here we have explicitly written out the positive, $Q Km and negative, 3>~(_ re ) m energy terms 
in the expansion of the eigenstates of H. The bag model states, $„ Km are eigenstates of 
K = /3(<r • L + 1) with eigenvalue 

^S m = -8gnP«C«- ( 7 ) 

The potential mixes P = + and P = — bag model states with opposite signs of k. 
In the bag model basis the Dirac equation is 

[(®«i>m\V\®LJ + (K* ~ sgnPE aK )8 p/>P 8 K > iK 8 al>a \ A p nKa = 0. (8) 

P',P,K',n,a',a 

The non-zero terms appearing in potential matrix element have (P', 
(+,+;+,+),(+,+;-,-),(+,-;+,-),(-,+;-,+),(+,-;-,+), and (-,-;-,-) and 
states obtained from these by (P ; , n') <->• (P, «). 

The P QK are the eigenstates of the free Dirac Hamiltonian E aK = p aK with momenta p aK 
fixed by the bag model boundary condition 

*t»(r)*i K (r)\ r =R = (9) 

where $ = $^7° and r is the bag radius. We consider only $+ K in Eq.Q since the negative 
energy states yield the same momenta p aK . 

We solve Eq.fjHJ) in the bag model basis up to some maximum value of the momentum 
Va m K- We include enough states in the basis so that the longitudinal sum rule, S^dql) (see 
Eq. of the next section) is satisfied. The value of the radius R is set by the requirement 
that radii of all included positive energy eigenstates of H have rms radii C E. R = 15 
fm is used to satisfy this condition. The negative energy eigenstates have a continuum of 
eigenvalues in the limit R — > 00 and therefore are P-dependent when R is finite. 

The diagonalization of Eq.fjHJ) is effected for each negative value of k. The spin symmetry 
admitted by H allows one to obtain the states with k > from those with k < since the 
states with k + k' = —1 are degenerate. For each value of k, angular integrations in the 
matrix elements of the potential in Eq.fjHJ) are evaluated explicitly. H is then diagonalized 



in a truncated basis of radial states (spherical Bessel functions appearing in the bag model 
states - see Appendix EJ) with wave numbers {p aK }^i- In this work we take a m = 600. 
Thus for each k we have 600 positive energy and 600 negative energy states in the basis. 
The eigenstates of Eq. (j3J) can be written 

where the angular part of the wave function is the spin-angle function, 3^ m , with orbital 
angular momentum I and spin-- coupled to total angular momentum j and projection m. 
The upper (lower) sign in the angular part of the lower component of the wave function is 
taken when k is negative (positive) with k = £ > and k = — (£ + 1) < 0. We obtain, from 
the expansion Eq.© and the explicit forms of the bag model wave functions Eqs. (|A2| ) and 
(|A3|) and their negative energy counterparts, the expressions for the radial wave functions 
(for k < 0) 



Mr) = [ A+ Ka N aK j e (p aK r) - iA n[ _ K)a N a( _ K) J - a{ k) — j e (p a{ _ K) r) ) (11) 



a=l 



9nn{r) = ( -^a^^ Y^-^je+i(PaKr) -iA n( _ K)a N a{ „ K) j £+1 (p a{ _ K) r)j. (12) 

Analogous expressions hold for k > 0. 

We note that the ground state has energy Eq = 840 MeV and rms radius (r 2 )^ = 0.66 
fm for string tension yfa = 1 GeV/fm. This string tension corresponds to a characteristic 
energy scale of cr 1//4 = 440 MeV. On the basis of dimensional analysis we find that Eq oc cj4. 



III. ELECTROMAGNETIC STRUCTURE FUNCTIONS 



The inclusive cross section for a spin-1 electromagnetic probe scattering from a hadronic 
target may be expressed to first order in the fine structure constant a as 
(fa 



dE'dQ 



a M (W 2 + 2VMan 2 (fl/2)) 



cm- 



Q 2 



W l + Wt [ l + 2^Uan 2 (#/2) 



Here <jm is the Mott cross section 
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a 2 cos 2 (0/2) 
AE 2 sin 4 (0/2) ! 
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FIG. 1: The longitudinal structure function W^"(|q|, v) vs. y = v — |q| for various values of the 
three-momentum transfer. 

with 9 the probe's scattering angle, E(E') its initial (final) energy and a q = Q g a, Q q is the 
charge on the hadronic target's constituent in units of the proton's charge and a is the fine 
structure constant. We have neglected the mass of the electromagnetic probe in the above. 
The inclusive longitudinal, Wl and transverse, Wt structure functions are related to the 
functions by 



r 



W 1 



W L = -Wi + 



Q 2 



Wo. 



(16) 
(17) 



|q| 2 — v 2 . 



We have assumed that g M = (z/, 0, 0, |q|) and defined Q = —q z 

The inclusive structure functions are determined from the matrix elements of the hadronic 
current on the target eigenstates and their energy eigenvalues as 



m,I 



w?(|q|^) = ^|(/]A + « + e^- r |o,m)| 2 5(£; 7 -£;o-^ 



(18) 
(19) 



m,I 



where |0, m) is the ground state with spin projection m = ±|, a + = 7°7+, 7+ = -^(7 1 +^7 2 ) ! 
and the \ J2 m i avera g es over initial spins and sums over final states, \I). We assume that 
the negative energy states are occupied and thus excluded from the sum over final states by 



Pauli exclusion. The + superscript on the structure functions indicates this exclusion. The 
projection operators on positive and negative energy states are defined by 

A + = E l J >< J l ( 20 ) 

I,Ej>E 

a- = £ ww* ( 21 ) 

I ) E I <E 

respectively. The states \I) are orthonormal (I\J) = Sjj. In the interest of subsequent 
discussion we also define W£ T with A_ in Eqs.flTHJ and (fTTJJ) instead of A + . 

The longitudinal and transverse structure functions W£ T are related to the response of 
the electromagnetic probe to the charge and spin fluctuations in the target, respectively. 
They are viewed as the response of the light, positive energy valence quark in a heavy-light 
meson, such as tu, as it is excited into a positive energy final state within the potential well. 
The longitudinal function W£ is associated with a spin-independent coupling of the probe 
to the target and thus the spin of the struck constituent is unchanged in contributions to 
it. The transverse function is associated with a spin-flip of the struck constituent and 
couples the upper and lower components of its wave function since 



7+ = _ , (22) 




where a + is the raising operator on Pauli spinors. 

We calculate the structure functions in a large basis of eigenstates of the Dirac Hamilto- 
nian [Eq.(jSJ] including values of k in the range — 100 < k < and < k < 100 in the sum 
over / in Eqs. (ll8| [TDJ) . We take 400 positive (negative) energy radial states each for each 
value of k for the calculation of W + {W~). We ensure that this basis supports all the states 
accessed by scattering up to |q| = 12.5 GeV to a level of 0.3% as measured by the static 
structure functions, S^dql) (as shown below). 

In order to obtain the structure functions as smooth curves we fold them with a Breit- 
Wigner function letting 

* (B/ - Bo _^IM__i__ (23) 

with level width T{y). We use the parameterization for T(v) given in Eq.(8) of Ref.Q with 
To = 100 MeV. The effects of hadronization of the excited states are taken into account 
qualitatively by this prescription. 
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Figure (JIJ shows (|q|, u) for various values of |q| plotted as a function of the variable 
y = v — |q|. We note some features of the W£. The curves at high |q| > 2.5 GeV, much 
larger than the characteristic energy scale of <r 1//4 = 440 MeV, exhibit scaling in the variable 
y. There is non-zero strength in the timelike region where y > corresponding to the fact 
that the confined constituent is bound in the target. The strength in the timelike region 
for |q| = 12.5 GeV is 7.3% and 11.5% for T = and 100 MeV, respectively. At |ql< 2.5 
GeV the curves exhibit resonances for y ~ — |q| and show Bloom-Gilman duality. 5j, |7| The 
quality of scaling is fair though scaling violations are significant in the timelike region having 
y > in the present |q| range. In order to understand these violations we will discuss sum 
rules. 

The full strength of the ground state is distributed among the the positive and negative 
energy eigenstates by the current 7 A V q ' r . Therefore both W£ and W£ must be taken into 
account to support the full strength of the ground state. If we write 

W L (\q\,v) = W L + (|qU) + W£(\ci\, u) (24) 

we obtain a function which satisfies the sum rule 

/oo 
dvW L (\<i\,v) = l, (25) 
-oo 

where S^Qql) is the longitudinal static structure function. Here we have integrated over all 
values of energy transfer v including the region v < and used completeness 

Ei j >< / i = 1 - ( 26 ) 

(In fact, in our numerical study, S^dql) has the value of unity to 0.3% or better for all 
considered |q|.) 

The total strength of the H^"(|q|, v) 

POO 

s£(|q|)= / w+(\ q \,v)<i (27) 

Jo 

since we have neglected the contribution of the negative energy states. As |q| increases, the 
strength S^dq)) due to the positive energy states decreases. It can be shown (see Appendix 
|B} that 

lim fi£(|q|) = i (28) 

|q|^oo L 
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FIG. 2: The transverse structure function WjT (|q|,y) vs. y = v — |q| for four values of |q| = 
2.5,5,10,12.5 GeV. 

At |q| = 12.5 GeV the S^Qql) = 0.52, independent of Tq. The strength missing from 
5^(|q|) is taken up by the negative energy states 

fi£(|q|)= f duWMv) (29) 

J — oo 

to maintain the sum rule Eq.(|25|). The scaling violations in W£ are due to the sharing of 
the strength of the ground state between the positive and negative energy eigenstates of the 
potential. 

We expect similar scaling violations in the case of the transverse structure function 
(|q|, v) which we have shown in Fig.(J2J). The strength in the timelike region for |q| = 12.5 
GeV is 1.4% and 4.5% for T = 0, 100 MeV, respectively. 

The transverse structure function Wr(|q|, u), defined analogously to Eq.pijl. satisfies the 
sum rule 

/oo 
duW T (\q\,u) (30) 
-oo 

= i^(0,m|(l-a z )|0,m) 

m 

= - + - = !, (31) 
6 6' K J 

where the first and second terms in Eq.()31|) correspond torn = — | and m = +|, respectively. 
The different strengths in m — ±| can be understood by analyzing the ground state wave 
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function. The matrix operator a + in Eq.(|19jl corresponds to a virtual photon with positive 
helicity and therefore can only scatter from components of the wave function which have 
Pauli-spin down. The ground state wave function is s-wave in the upper component and 
the photon can't scatter from it when m — +|. It only scatters from the spin down part 
of the p-wave lower component resulting in less strength for the m — +| state than in the 
m = — \ state. 

In deriving the result Eq. fpITj) we used the facts that 

dr r 2 / 2 (r) = \ (32) 
dr r 2 g 2 (r) = ± (33) 

where r/o(r) and rgo(r) are the upper and lower radial wave functions for the ground state. 
These follow directly from the Dirac equation for a massless constituent independent of the 
string tension \fa as illustrated in Appendix O Additionally it can be shown that 

E = 2(V) , (34) 

where (O) = J d 3 r^l(r)O^ (r). 

Up to this point we have calculated the structure functions W£ T for a confined Dirac 
valence particle. These functions concern the scattering of the valence particle in the ground 
state to a positive energy excited state. The contribution of the negative energy states in the 
sum over intermediate states I in Eq. (jl8|) has been ignored since these states are assumed 
to be filled. 

We now consider the response of the vacuum in the presence of the infinitely massive 
color source. The response of the vacuum corresponds to scattering from a negative energy 
state into a positive energy state. Consider the contribution to the vacuum response when 
the particle is initially in a negative energy state / and scatters into the unoccupied ground 
state: 

^| 0) ~(|qU) = ^|(0,m|e-^ r A_|/)| 2 5( J E; / - J E;o + ^). (35) 

m,I 

Here the superscript (0)— indicates that we are scattering from a negative energy state 
into the ground state. We have taken the energy transfer v —>■ — v since a negative energy 
state is equivalent to a positive energy antiparticle. This contribution is shown Fig.® as a 
function of y for |q| = 10 GeV. We noted above that the negative energy states depend on 
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FIG. 3: The vacuum contribution to the longitudinal structure function (|q| , v) for scattering 

from a negative energy state into the ground state vs. y = u — |q| for |q| = 10 GeV. Only the solid 
curve in the region y < is physical. The continuation (dashed curve) into the region y > has 
no meaning. 

the cavity radius, R, since they are not bound by the potential Eq.fjHJ). However is 
only weakly dependent on the cavity radius since the matrix element receives contributions 
for r < 1 fm, the radial extent of the ground state. At these radii the negative energy 
states are not sensitive to the boundary condition. This was verified for R — 10 and 15 
fm. We interpret as the contribution to the vacuum response in the spacelike y < 

region. It is a consequence of the fact that interactions move some of the strength of the 
vacuum response from the timelike region into the spacelike region. We cannot reliably 
calculate the response in the timelike region since confinement acts at all distances in our 
model. Additional contributions to the vacuum response in the spacelike region correspond 
to scattering from a negative energy state into a positive energy state J and are denoted 
. These contributions are smaller than the since the first excited state has 

energy Ej = i = 1100 MeV, 260 MeV above the ground state. 

IV. PLANE WAVE IMPULSE APPROXIMATION 

We may compare the exact structure functions calculated in the previous section with 
those obtained in PWIA. We consider only the longitudinal structure function in this section. 
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FIG. 4: The longitudinal structure function W£"(|q|, i/) vs. y = v — |q| in PWIA and the exact 
result for |q| = 10 GeV. 

In PWIA, one assumes that the state of the struck constituent may be described as a plane 
wave with energy 

£ k+q = |k + q| + <y) , (36) 

for a massless constituent. The expectation value of the potential [Eq.®] in the ground 
state, (V)q is required to give the correct value for the energy weighted sum rule 

POO 

4 1)+ (|q|)= / dvvW+{H,v). 



(37) 

We may calculate this quantity analytically in the limit |q| —>■ oo using the technique in 
Appendix IH1 to obtain 



lim 4 1)+ (lql) , 

|q|^oo Z 



q 



(38) 



The PWIA neglects interactions of the struck constituent in the final state. The expres- 
sion Eq.(|18|). for example, then simplifies to 



W 



L,PWIA 



d 3 k , ,9 

(M k+q ,,|e iq - r |0, m) <J(|k + q| + (V> -Eq-u), (39) 



(27T) 



where |uk+ q , s ) is the positive energy free-particle Dirac spinor with spin s. An analogous 
expression holds for the W^ PWIA . 

Evaluating the matrix elements in the above equation gives 



W 



If d 3 k 
L.rwiAW ~ 2 7 4tT 

x 5(k z -y+ (V) - E ) 



f^k) + g' (k) - 2f (k)~g (k) 



(40) 
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where we have neglected terms C(j^j) and fo(k), go(k) are the Fourier-Bessel transforms of 
the radial wave functions appearing in the Fourier transform of the ground state: 

3 / fo(k)y°i (k) \ 
*o,«(k) = (2*)i *> m ): . (41) 

The y l j m are the spin-angle functions. Figure 0] shows the PWIA approximation for the 
longitudinal structure function and the exact result. 

The argument of the 5 function depends only on the variable y and therefore the PWIA 
response also exhibits scaling behavior. It has scaling violations which are ^(■yO f° r IqI 
finite. The W£ PWIA satisfies the sum rules Eq.(J2BJ) and Eq. lfHTJ) but has a different shape 
than the exact curve. It has more response in the timelike region than the exact result 
with 14.7% of the total strength in the timelike region compared with 11.5% for the exact 
result. This is a consequence of the interference of the upper and lower components of the 
ground state wave function in the third term in square brackets of Eq. (|40j) . Neglecting this 
term, the W^ PWIA would peak at y = — nearly coincident with the peak of the exact 
curve. The interference term shifts strength into the timelike region. In a subsequent work 
we shall explore the consequences of such effects on spin-dependent observables such as the 
helicity structure function, jlj We note finally that the PWIA and exact curves in Fig.(j3J) have 
different widths. The larger width of the exact curve is due to FSI of the struck constituent 
with the potential, V. As in the case of the semirelativistic model |5], FSI effects persist 
in the limit |q| — > oo. We note that Brodsky and collaborators have shown that FSI affect 
the interpretation of parton distribution functions [^| and can contribute to large single spin 
asymmetries in semi-inclusive DIS processes jlfll ]. 



V. CONCLUSION 



We have calculated exactly the eigenstates of a massless Dirac particle confined in a linear 
potential for excitation energies up to ~ 10 GeV. The unpolarized valence quark structure 
functions of the particle in the ground state scattering into a positive energy excited state 
were determined by coupling a spin-- electromagnetic probe to the conserved vector current 
of particles in the potential well. The structure functions exhibit y or Nachtmann scaling, 
a generalization of Bjorken scaling. The exact response takes account of interactions of 
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FIG. 5: The deviation from Callan-Gross relation in terms of the function 5 of Ea. (|42j) for |q| = 10 
GeV. The dashed curve is the scaled model calculation (see text). The solid curve is the MRST2002 
pdf fit to data. 

the confined particle both before and after interaction with the electromagnetic probe and 
therefore y or Bjorken scaling is not a consequence of treating the constituent of the target as 
a free particle. In contrast to Refs. [ill . Il^ who find that FSI have no effect on the structure 
functions, and therefore also conclude that Bjorken scaling applies to systems of interacting 
constituents, the present model does not reproduce the parton model results when FSI are 
taken into account. 

Another example of the role interactions have in determining the structure function of 
confined Dirac particles can be seen if we consider the function 



S(y, |q|) = 



F 2 (y, 


q\)-2xF 1 (y, 


q|) 


F 2 (y, 


q\) + 2xF 1 (y, 


q|) 



(42) 



where F\ = MW\ and F 2 = z^W^, M is the mass of the target. At leading order in the 
quark-parton model one obtains 5(y, |q|) = (the Callan-Gross relation [l^l) for all y and 
|q|. In Fig.© we have plotted, for |q| = 10 GeV, the MRST2002 |lj| fit (solid curve) to 
deep inelastic data using the next-to leading order prediction jlo| 



where < x < 0.4 and g(x, Q 2 ) is the prediction for the gluon. It is compared to the present 
model calculation (dashed curve) which has been rescaled in order to obtain agreement with 




(43) 
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the first moment of the up valence quark distribution for the MRST pdf assuming that the 
response of the proton is due to a single up valence quark. The non-perturbative effects of 
interactions could be the origin of the difference in the two curves in Fig.flSJ). 

In future work we will apply the present model to the study of spin-dependent effects 
such as the helicity and transversity distribution of the constituents of the target. We expect 
that these observables depend strongly upon the relativistic treatment of the constituent as 
a Dirac particle. The comparison of the present model with experimental data requires one 
to include radiative gluon effects to obtain the Q 2 evolution of the structure functions. This 
problem is currently under study. 



APPENDIX A: BAG MODEL BASIS STATES 



The free Dirac Hamiltonian Hq is 

[a • p + pm] § p aKm = sgnP^jw) 2 + m 2 <$> p aKm 

= sgnPE aK ^ Km , (Al) 

for a particle of mass m (not to be confused with the total spin projection which appears as 
the subscript m in the wave function $„ Km ). The complete set of states includes states of 
positive E aK > energy states, 'particles' and negative E aK < energy states, 'antiparticles.' 
The bag model states are solutions of Eq. (jAl|) within a cavity of radius R. These particle 
states may be written as 



K Km (r) =N a J . r-^ aK , , k > 0. (A3) 



Here the normalization N aK is fixed by the condition f d 3 r$+l m (r)$+ Km (r) = 1 to give 



N aK = { R 



2(jE ^ )2 + —L-( m + 2E aK K) 



j±(«+i)(P^) \ ■ (A4) 



. (Pa K ) 2 (Pa K ) 2 R 

The upper (lower) sign corresponds to k positive (negative). 

The antiparticle (negative energy) states are obtained from Eqs. f)A2|) and (jA3|) by in- 
terchanging upper and lower components and multiplying one of them by —1. With this 
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prescription the sign of k for the negative energy state is opposite to the sign of k of the 
positive energy state from which it was obtained. 

As mentioned in the text, we determine the allowed p Ka by subjecting the $ QKm to the 
bag model boundary condition at r = R 

$®\ r=R = 0, (A5) 

where $ = &(3. Substitution of the wave functions in Eqs. (jA2| ) and (jA3|) into this expression 
results in the following transcendental equations; for k = — (£ + 1) < 



Je(p aK R) = \/ ^ aK , m je+i(PaKR), (A6) 
E aK + m 



and for k = I > 



3t{p aK R) =-\ , m je~i(PaKR)- (A7) 

V E aK + m 

The subscript a = l,...,oo indexes the solutions of the transcendental equations. The 
number of nodes in the upper component radial function is (a — 1). 



APPENDIX B: Sj(|q|) IN THE LIMIT |q| -» oo 

The expression for the ^(Iql) is given by 

poo 

5+(|q|) = / <W+(|q|) (Bl) 
Jo 

= ^^(0,m|e^ q - r A + e iq - r |0,m). (B2) 

m 

We may evaluate this expression analytically in the limit |q| — > oo by recognizing that in 
this limit the intermediate states \I) in Eq. (|B2|) may be taken to be free-particle solutions, 
|^k+ q ,s) of the Dirac equation since their overlap with the ground state is dominated by 
small radii < 1 fm where the exact wave functions resemble free-particle wave functions. 
Thus we replace 

e -«,r A+e «, r ^ J _0_| Mk+q)S )( Mk+q;S | e ^ r , (B3) 



valid in the limit |q| — » oo and obtain upon substitution into Eq. (jB2 



1 f d k 
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The projection operator for free-particle solutions to the Dirac equation, expressed to 

°(w)- is 

Et 1 ( a ' k± \ 

Wk+ q ,s< +q , s = - [l + a z + — — , (B5) 



2 V |q| 

where q = |q|z and kj_ = (k x ,k y ). Substitution of Eq. (jB5|) into Eq. ()B4|) gives the result 
Eq.fEHD- 

APPENDIX C: GROUND STATE NORMALIZATIONS 

The Dirac equation for the ground state gives 

f (r) = E g (r) (CI) 
-g'oir) - \{r) = (E - v^r)/o(r). (C2) 

for the upper and lower radial wave functions. Multiplying Eq. (|ClJ) by go(r) and Eq. fj( J2j) 
by fo(r), summing and integrating gives 

poo i*oo 

E / drr 2 {f 2 Q {r) - g 2 (r)) = ^ \ drr 2 f*(r)r. (C3) 
Jo Jo 

The r.h.s. of the above is the expectation value of the potential V(r) = ^/ar{l + (3) in the 
ground state since it doesn't couple to the lower components of the wave function. Substitute 
Eq.flHU) into Eq.JCSl) to obtain 

J Q drr\f (r) - g*(r)) = -. (C4) 

The normalization condition is: 

POO 

/ drr\ft{r)+gl(r)) = l. (C5) 

Adding and subtracting these equations gives the result Eqs.(j32J), independent of the string 
tension, ^/a. 
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